Theory of Spin Relaxation in Two-Electron Lateral Coupled Quantum Dots 
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A global quantitative picture of the phonon-induced two-electron spin relaxation in GaAs double 
quantum dots is presented using highly accurate numerics. Wide regimes of interdot coupling, mag- 
netic field magnitude and orientation, and detuning are explored in the presence of a nuclear bath. 
Most important, the giant magnetic anisotropy of the singlet-triplet relaxation can be controlled 
by detuning switching the principal anisotropy axes: a protected state becomes unprotected upon 
detuning, and vice versa. It is also established that nuclear spins can dominate spin relaxation for 
unpolarized triplets even at high magnetic fields, contrary to common belief. 



Electron spins in quantum dots [ij are among perspec- 
tive candidates for a controllable quantum coherent sys- 
tem in spintronics 0, Q ■ Spin qubits in GaAs quantum 
dots, the current state of the art 0, [H, are coupled to 
two main environment baths: nuclear spins, and phonons 
0. The nuclei dominate decoherence, which is on fis 
timescales. But only phonons are an efficient energy sink 
for the relaxation of the energy resolved spin states, lead- 
ing to spin lifetimes as long as seconds [7j. 

The extraordinary low relaxation is boosted by or- 
ders of magnitude at spectral crossings, unless special 
conditions — such geometries we call easy passages — are 
met [1, 0| • Spectral crossings seem inevitable in the ma- 
nipulation based on the Pauli spin blocka de [H lioj , the 
current choice in spin qubit experiments On the 

other hand, a fast spin relaxation channel may be de- 
sired, e.g., in the dynamical nuclear polarization 12h14|. 

The single-electron spin relaxation is well understood 
[l5l[l6]: it proceeds through acoustic phonons, in propor- 
tion to their density of states, which increases with the 
transferred energy. The matrix element of the phoiron 
electric field between spin opposite states is nonzero due 
to spin-orbit coupling or nuclear spins. At anticrossings 
the matrix element is enhanced by orders of magnitude, 
even though the anticrossing gap is minute /^eV) . The 
relaxation rate can be either enhanced or suppressed, de- 
pending on whether the energy or the matrix element 
effects dominate. 

The two electron relaxation rates were measured in 
single [iMl and in double dots 



works so far mostly focused on single dots [23|, l24 | 
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vertical double dots |25l . |26[ , in which the symmetry of 
the confinement potential lowers the numerical demands. 
A slightly deformed dot was considered iir Refs. [27, 28], 
and a lateral coupled double dot in silicon in Ref. [29|. 
What is key for spin qubit manipulation and most rele- 
vant for ongoing experiments, is the case of weakly cou- 
pled and biased coupled dots. In addition, the relative 
roles of the spin-orbit and hyperfine interactions in the 
spin relaxation in GaAs quantum dots has not yet been 
established. 

The analysis of the two electron double dot relaxation 



is challengiirg because many parameters need to be con- 
sidered simultaneously: the magnitude and orientation of 
the magnetic field, the orientation of the dot with respect 
to the crystallographic axes, the strength of the interdot 
coupling (parametrized by either tunneling or exchange 
energy) and the bias applied across the double dot (de- 
tuning). Here we cover all these parameters, including 
the nuclear bath, providing specific relevant predictions 
for experimental setups [30|. Perhaps the most striking 
results are the existence of islands of iirhibited spin re- 
laxation in the magnetic field and detuning maps, and 
the switch of the two priircipal C2v axes along which the 
relaxation shows a minimum or maximum, as detuning 
is turned on. While singlets and polarized triplets relax 
by spin-orbit coupling, the spin-unpolarized triplet relax- 
ation is dominated by nuclear spins over a wide param- 
eter range (the spin-orbit induced anisotropy is wiped 
out), contrary to common belief. The predicted giant 
spin relaxation anisotropy is a unique and experimentally 
testable signature of spin-orbit spin relaxation, which can 
also be useful for spin nanodevices, as we argue in this 
paper. 



MODEL 

We consider a laterally coupled, top-gated GaAs dou- 
ble quantum dot patterned in the plane perpendicular to 
z = [001]. In the two-dimensional and envelope function 
approximation, the Hamiltonian reads 



H = J2 {T, + V, + Hz,, + H,. 

i=l,2 



^)+Hc, (1) 



where i labels electrons. The single-electron terms are 



T = PV2m = {~ihV + eAy /2m, 



(2) 



V = (1/2) mwgmin{(r - d)^ (r -I- d)^} eE • r, (3) 



Hz = (,g/2)^BT-B, 

Hso = Hhi + Hd + Hd3, 



(4) 
(5) 



the kinetic energy, the biquadratic confinement poten- 
tial, the Zeeman term, and the spin-orbit couplings, re- 
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spectively. The position and momentum vectors are two- 
dimensional, where x = [100] and y = [010]. The proton 
charge is e and the effective electron mass is m. The con- 
finement energy, Eq = Hluq, and the confinement length, 

1/2 

^0 = {h/moJo) , define the characteristic scales. The 
potential is minimal at ±d and we call 2d/lo the interdot 
distance. The electric field E is applied along the dot 
main axis d. Turning on E shifts the potential minima 
relative to each other by the detuning energy e = 2eEd. 
The magnetic field is B = {Bx, By,Bz). We use the sym- 
metric gauge, A = Bz{-y,x) /2, and a = {(Jx,Uy,(Jz) 
are the Pauli matrices. The Lande factor is g and the 
Bohr magneton is /ie- The Bychkov-Rashba, and the 
linear and cubic Dresselhaus Hamiltonian read 

Hh, ^ (?l/2mZbr) {(^xPy - CTyPx) , (6) 

Hd = {h/2mU) {-a^P^ + ayPy) , (7) 

i?d3 = {jc/2h^) {cTxPxPy - TyPyPl) + H.C., (8) 

parameterized by the spin-orbit lengths Zbr and /d, and a 
bulk parameter 7^. Nuclei, labeled by n, couple through 



(9) 



where /3 is a constant, and I„ is the spin of a nu- 
cleus at the position R„. The Coulomb interaction is 
He = e^/47re |ri — r2|, with the dielectric constant e. 
The Hamiltonian, Eg. (U), and its energy spectrum is dis- 
cussed in Refs. |3ll |32|. including our numerical method 
(configuration interaction) for its diagonalization. Here 
we extent it by including nuclear spins, which we treat by 
averaging over unpolarized random ensemble. See Sup- 
plementary [33| for further details. 

The relaxation is mediated by acoustic phonons 
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FIG. 1. Calculated energies of the lowest states for (a) vari- 
able interdot coupling (at B = 5 T), and (b) detuning (at 
B = 2T). Singlet states are given by dashed, triplets by solid 
lines. The blue strokes mark singlet-triplet anticrossings. In 
(a), the energy of To is subtracted, and in (b), the quadratic 
trend in E is subtracted. The green arrows denote points 
of exact compensation and the red oval in (b) shows where 
nuclear spins dominate the To relaxation. 



where = (i|e'*^'^|j) is the matrix element of the 
states with energy difference hiOij. Here we are inter- 
ested in the rates of the singlet {S) and the three triplets 
(T-|_, Tq,T^) at the bottom of the energy spectrum. 

In numerics we use GaAs parameters: m — G.067TOe, 
with me the free electron mass, g = —0.44, q = 5290 
m/s, ct = 2480 m/s, p = 5300 kg/m^, de = 7 eV, ehu = 
1.4 X 10^ cV/m, e = 12.9, 7^ = 27.5 eVA^, /3 = 2 fieY 
nm^, 1=3/2. We choose typical lateral dots values, ^br = 
2.42 /im, ~ 0.63 /um, d || [110] and the confinement 
energy £^0 = 1-0 meV, corresponding to Iq = 34 nm. 
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(10) 

= cTe, and piezoelectric potentials, 

i'iiqxqye^^, + qzqxeQ^y+qyqze^^;,)/Q^- The 
phonon wave vector is Q, and the electron position vec- 
tor is R = (r, 2). The polarizations are given by A, the 
polarization unit vector reads e, and the phonon annihi- 
lation (creation) operator is denoted by b (6^). The mass 
density, the volume of the crystal, and the sound veloci- 
ties are given by p, V, and c\, respectively. The phonon 
potentials are parameterized by Ce, and /114. 

We define the relaxation rate as the sum of the in- 
dividual transition rates to all lower-lying states for 
both piezoelectric and deformation potentials. Each rate 
(from |i) to |j)) is evaluated using Fermi's Golden Rule 
in the zero-temperature limit. 



(11) 



We start with an unbiased double dot. We plot its 
spectrum in Fig. [1^) as a function of the interdot cou- 
pling, which translates into an exponentially sensitive 
5 — To exchange splitting J. Electrical control over J, 
necessary e.g. to induce the vSWAP gate ^Ij], allows for 
a fast switching between the strong and weak coupling 
regime, corresponding to the exchange splitting being 
larger and smaller than the Zeeman energy, respectively. 
During this switching, the ground state changes at an 
S — T+ anticrossing. 

We cover the freedom of the interdot coupling in Fig. [2] 
Panel a) shows the relaxation of the first excited state [S 
or r+, see Fig. [1^)]. First to note is the strong relaxation 
suppression at the S — T+ anticrossing as the transferred 
energy becomes very small. Remarkably, the anticross- 
ing does not influence the rate of Tq, plotted at panel b), 
at all (the peak close to d = is due to an anticrossing 
with a higher excited state). Even though the dominant 
channel, Tq ^ is strongly suppressed here, its reduc- 
tion is exactly compensated by the elsewhere negligible 
Tq ^ S channel. The exact compensation arises for the 
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FIG. 2. Calculated relaxation rates of (a) the first excited 
state {S or T+, see Fig.[T^)) and (b) the triplet To as a function 
of the in-plane magnetic field orientation 7 — a,Tccos{Bx / B) 
(angle) and the interdot distance 2d/lo (radius of the polar 
plot), for a double dot at _B = 5 T. The x and y axes cor- 
respond to crystallographic axes [100] and [010], respectively. 
The dot orientation d || [110] is marked by a line. The blue 
half circles indicate the S — T+ anticrossing, also marked on 
Fig. [1^). The x axis is converted to the tunneling energy T 
and the exchange J, in addition to 2d/lo. The rate is given 
in inverse seconds by the color scale. The system obeys C2v 
symmetry, so point reflection would complete the graphs. 
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FIG. 3. Calculated relaxation rates of (a) the first excited 
state and (b) the triplet To as a function of the in-plane mag- 
netic field orientation 7 (angle) and the magnetic field magni- 
tude (radius of the polar plot), for a double dot with T = 0.1 
meV. The layout with respect to the crystallographic axes is 
the same as in Fig. 2. The rate is given in inverse seconds by 
the color scale. 



the coordmates along the dot axes Xd,yd = {x ± y)/V^, 



relaxation into a quasi-degenerate subspace (we denote 
such cases on Fig. 1 by green arrows) if 

AE <$:mm{E,hcx/lo}- (12) 

Here E is the transition energy and AE is the energy 
width of subspace (the anticrossing gap) . Equation 
states that the energy width AE is too small to be re- 
solved by either phonons with energy E or electron wave 
function scale Iq (33| . The relaxation then proceeds into 
the subspace rather then into its constituent states, so 
that any mixing of the states within the subspace is ir- 
relevant. 

Further to note on Fig. 2 is the anisotropy of re- 
laxation, which reflects the anisotropy of the spin-orbit 
fields. In the weak coupling regime, the relaxation rates 
are minimal if the magnetic field orientation is parallel to 
the dot main axis, which results in an isle of strongly pro- 
longed spin lifetimes. Note that this is in contrast to the 
biased dot (see below), and to the single-electron case, 
where the minimal in-plane magnetic field direction, the 
easy passage, of a d || [110] double dot is perpendicular 
to d [3, 1131 • The switch can be understood from the efi^ec- 
tive, spin-orbit induced, magnetic field 9;] if written using 



B,o - B X {xdii^; -i^')[no] + y4i-; + i-')[no]}/V2. 

(13) 

At the anticrossing, the mixing due to is by far domi- 
nant, so the minimum appears with B along [110]. This 
Xd dominance will be the case for a biased dot, too. On 
the other hand, in a single dot Xd and yd induce com- 
parable mixing, and Bso becomes minimal if the larger 
term (the one with yd) is eliminated. Weakly coupled 
unbiased dot is in this respect similar to a single dot as 
the two-electron transitions can be understood as flips of 
a particular electron located in a single dot. Since the 
direction for the rate minimum switches upon changing 
d, the system does not show an easy passage, that is a 
low relaxation rate from weak to strong coupling regime. 

We plot the magnetic field dependence for a weakly 
coupled unbiased double dot in Fig. |3] and observe similar 
behavior as in Fig. 2. The relaxation rate is minimal if 
B 1 1 d throughout the shown parametric region. This 
is because the anticrossing and the related directional 
switch happens here at so small magnetic field that it is 
not visible at the figure resolution. For completeness, we 
note that the r_ relaxation behavior is very similar to 
the one for Tq on both Figs. 2 and 3, and we do not show 
it. 
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BIASED DOUBLE DOT 

We now consider a biased double dot. Its spectrum is 
shown in Fig. [TJd) as a function of the detuning. The 
ground state singlet is in the (1,1) configuration (one 
electron in each dot) for low, and in the (0,2) configu- 
ration (both electrons in one dot) for large detunings. 
The crossover, a broad singlet-singlet anticrossing, is a 
key handle in spin measurement and manipulation [ll| . 
The low to large detuning crossover involves S — Tj- an- 
ticrossing, exploited for nuclear-spin pumping 

We show the detuning and magnetic field influence on 
the relaxation in Fig. |4l At the singlet-triplet anticross- 
ings, we observe that first, the relaxation rate of the first 
excited state dips at the S — T+ anticrossing (though the 
dip is very narrow and hard to see at the figure reso- 
lution), and second, the T_ rate strongly peaks at the 
S — T_ anticrossing. This is a demonstration of the dom- 
inant effect of the anticrossing on the transition energy, 
and matrix element, respectively. Third, there are no 
other manifestations of the S — T± anticrossings, a fact 
due to the exact compensation already mentioned before. 
The anisotropy features of this geometry are striking. In 
the given range of detuning energies, states except Tq ex- 
hibit a very distinctive easy passage for a magnetic field 
along [110], where the relaxation is up to to three orders 
of magnitude smaller than with B along [110]. Though 
the directional switch occurs — rates become minimal for 
a magnetic field along [110], it is again out of the fig- 
ure scope (very small and very large detunings). The 
rates increase at detunings > 2 meV, because of spectral 
crossings with excited triplets. Fig. [T]d), regime normally 
avoided in experiments. Double dots, with their spectral 
idiosyncrasies, are a unique system to observe a giant 
amplification of the spin-orbit anisotropics by a physical 
observable with bias control. 

In large parts of the parametric space the relaxation 
of To is dominated by nuclear spins, thus being isotropic. 
This is surprising, since the effective (Overhauser) nu- 
clear magnetic field Bnuc is of the order of mT, much 
smaller than the spin-orbit field in Eq. ([13]), Bgo ^ 
{lo/lso)B ~ 30 mT at i? = 1 T for our parameters. One 
therefore expects the nuclei to lead to much slower re- 
laxation than the spin-orbit coupling. This was indeed 
the case for the unbiased dots and Figs. 2 and 3. How 
can then nuclei dominate here? Looking on Fig. [T|d, this 
happens when states Tq and ^(l, 1) are nearby in energy. 
Here, the otherwise negligible hyperfine effects take over, 
because the s pin - orbit induced mixing of these two states 
is forbidden 27[. Estimating the wave function admix- 
ture in the lowest order, the nuclei dominate if 



Bso/\Eto — Ek\ < Bnuc/\ETo — Es\ 



(14) 



with k being the closest state to which Tq is coupled by 
the spin-orbit interaction. The above condition general- 
izes in an obvious way for other states than Tq and there 
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FIG. 4. Calculated relaxation rates of (a) the first excited 
state, (b) To, and (c) T_ as a function of the in-plane magnetic 
field orientation 7 (angle) and detuning energy (radius of the 
polar plot), for a double dot with 2d/lo = 4.35 (T = 10 fieV), 
chosen along Ref. and B = 2 T. The layout with respect 
to the crystallographic axes is the same as in Fig. 2. The 
rate is given in inverse seconds by the color scale. The blue 
lines indicate the singlet-triplet anticrossings, in line with the 
marks in Fig. [TJj). The dashed red lines in panel b) confine 
the area where hyperfine coupling dominates. 



are additional cases of nuclear dominance in our system. 
However, they happen on parameter regions too small to 
be visible on the resolution of Fig. 4, so we discuss them 
only in the Supplementary material [33| . 



CONCLUSIONS 

Our predictions are experimentally observable. Until 
now the spin-orbit origin, and especially its induced di- 
rectional anisotropy of the spin relaxation in weakly cou- 
pled two-electron dots has not yet been experimentally 
established. With employing vector magnets it should 
now be possible to overcome earlier experimental chal- 
lenges and change the magnetic field orientation while 
keeping the sample fixed and detect the anisotropy [36| . 
The spin-orbit / nuclear induced relaxation can be masked 
by cotunneling and smeared by a finite temperature. The 
former is reduced in the charge sensing readout setups 
37|, in which the coupling to the leads can be made 
small. The latter effect is small for experimentally rele- 
vant sub Kelvin temperatures, such that the directional 
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anisotropies are well preserved. 

Our results demonstrate control over the spin-orbit 
induced anticrossing gaps (easy passages appear if the 
gaps are closed) by sample and magnetic field geome- 
try. It offers electrical tunability of spin relaxation, by 
changing the double dot orientation (in the Supplemen- 
tary material [s^, we suggest a spin current measure- 
ment device exploiting easy passage). In addition, such 
control may be especially useful when dealing with hy- 
perfine spins. Indeed, in the polarization scheme con- 
sidered in Ref. the nuclear spin polarization is pro- 
portional to non-hyperfine assisted spin relaxation (see 
Eq. (7) therein) and so would benefit from a setup with 
maximized spin-orbit induced relaxation rates (out of the 
easy passage). On the other hand, the adiabatic pumping 
scheme demonstrated in Ref. [s^ , relies on the an- 
ticrossing being solely due to the nuclear spins (and not 
the spin-orbit coupling), suggesting improved efficiency 
in an easy passage configuration. We propose a similar 
non-adiabatic nuclear pumping scheme based on the easy 
passage in the Supplementary material [11] . All these ex- 
amples illustrate the potential benefits which intentional 
control of spin relaxation, based on our results, may offer. 
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SUPPLEMENTAL INFORMATION 

In this supplementary material, we provide further 
details illustrating the main text, and derivations of 
some of its results. In Sec. A we comment on our nu- 
merical method and derive the discretized form of the 
electron- nuclear Hamiltonian, Eq. (9). In Sec. B we de- 
rive Eq. (12), the condition for the exact compensation 
and illustrate the exact compensation showing channel 
resolved relaxation rates. In Sec. C we compare the nu- 
clear vs. spin-orbit induced relaxation rates and discuss 
the additional cases to the one mentioned in the main 
text, where the nuclear spins dominate the relaxation. 
Finally, in Sec. D we suggest schemes for dynamic nuclear 
polarization and detection of spin polarization, which are 
based on easy passages. 

Numerical method 

We use the exact diagonalization of Eq. (1) in the con- 
figuration interaction method. In this work, the two- 
electron basis consists of 1156 Slater determinants, gen- 



erated by 34 single electron orbital states. The discretiza- 
tion grid is typically 135 x 135. The relative error for en- 
ergies is below 10~^. The diagonalization procedure was 
described in detail in Ref. [32] . Here we extend it adding 
nuclear spins, for which we now derive the discretized 
form of the Hamiltonian in Eq. (9), H^^^. 

Consider a basic element of the spatial grid, a rectan- 
gular box with lateral dimensions hx and hy. Such vol- 
ume elements are labeled by the index k — 1, M , with 
M their amount. In the two dimensional approximation 
one assumes that the electron wave function along the z 
direction is fixed to ip{z). Let it be, for concreteness, the 
ground state of a hard- wall confinement of width w — 11 
nm. As follows from Eq. (9), the discretized Hamiltonian 
is diagonal in the spatial index. It has matrix elements 

{H^ll^)kk' = Skk'^Y^^n ■ fT{k,^\S{r - 'R„)\k,i^), (15) 

where the sum is over nuclei inside the volume element 
k. To proceed further, we discretize the delta function in 
two dimensions as 1/hxhy, introduce the nuclei volume 
density as I/vq, and get 

{H^nc)kk' ^ 6kk'{P/vo)Y,voln- n\^}:{z^)\yKhy. (16) 

We now replace the sum over (typically many) nuclear 
spins by an effective spin X and get the discretized form 
of the Hamiltonian as 

{H^:^)kk'=5kk'{P/vo)I-<r. (17) 

By the central limit theorem, the effective spins are com- 
pletely described by their average and dispersion, which 
follow from the corresponding characteristics of the nu- 
clear spin ensemble. For random unpolarized nuclear 
spins, which we consider, it holds 

(I„)=0, (I„-I„) =5„„J(/+1), (18) 

so that the effective spins have zero average and the fol- 
lowing dispersion 

{Ik-Ik')=&kk'I{I + l)/N. (19) 

Here N is the number of nuclei in the grid volume el- 
ement, N = hxhyhz/vo, where the effective extension 
along z is defined by the wave function profile 

h-' - I dzmz)\\ (20) 

For the hard- wall potential one gets hz = 2w/3. 

To calculate relaxation rates, we average Eq. (11) over 
typically 50 configurations of nuclear ensemble. A sin- 
gle such configuration is parametrized by a set of effec- 
tive spins X drawn from a random Gaussian ensemble 
described by Eq. (A5), for which we diagonalize the two- 
electron Hamiltonian. Having the two-electron spectrum 
allows us to calculate the matrix elements Mij and en- 
ergy differences ujij which enter the relaxation rates Fy 
in Eq. (11). 
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Exact compensation of relaxation rate channels 

In general, the relaxation rate channels significantly 
change at spectral anticrossings because of the strong 
mixing of states. We consider the total relaxation rate 
by summing over the individual relaxation channels of all 
lower lying states. Therefore, a change in one relaxation 
channel may be compensated by another channel, such 
that the total relaxation rate is smooth (no peak or dip) 
across the anticrossing. This generally happens if a state 
relaxes into a quasi-degenerate subspace of anticrossing 
states. We exemplify such a situation on Fig. lb) by a 
green arrow and consider relaxation of Tq (state i) at this 
detuning. We write the total relaxation of Tq as 

where the sum includes the two quasi-degenerate states 
ji = T_|_ and j2 = S. The exact compensation arises if 
one can approximate the argument of the delta function 
by a common energy difference oj^^ ujij^ uj. Indeed, 
the relaxation can be then written as 

r» = ^ E ^ I^Q.aI' mPMW{^ - ^q), (22) 

^ Q.A ^ 

where Pj = is the projector on the quasi- 

degenerate subspace, which is not influenced by mixing 
of the basis states j. The condition for the approxima- 
tion Eq. (B2) to be valid is that both under-integral fac- 
tors in Eq. (Bl), the phonon density of states as well as 
electron overlap integral, are not changed much by the 
slight shift of the transition energy. The phonon density 
of states scales as a certain power (albeit different for 
piezoelectric and deformation potential) in the phonon 
wave vector, which translates into the condition 

Wjija/CA < Wij/CA. (23) 

On the other hand, the natural scale for the electron wave 
function is the confinement length ^o, so that the overlap 
M will not change much if 

^nj2/cx « l/lo- (24) 

Denoting ujjij2 = AE/h we get Eq. (12) of the main 
text. For the interdot coupling denoted on Fig. la and 
the dominant piezoelectric phonons we have E = 125 
/xeV, hct/lo = 48/zeV, while AE is just 7/zeV, so that 
the exact compensation condition is well satisfied. 

To illustrate the exact compensation, we plot in Fig. [5] 
the individual relaxation channels as a function of inter- 
dot distance. The parameters are chosen the same as in 
Fig. 2 of the main text. We find the exact compensation 
at the S — anticrossing for the Tq and the T_ relax- 
ation. In the case of the unpolarized triplet, the dip of 




interdot distance Sd/l^ 

FIG. 5. Calculated channel resolved relaxation rates vs. in- 
terdot distance in units of Iq for both parallel (top) and per- 
pendicular to d (bottom) in-plane magnetic field orientation 
(B = 5 T, zero detuning). The relaxation channels of To and 
T- are in blue and black color, respectively. The relaxation 
rate of the first excited state is red. 

the Tq — )■ r+ channel is compensated by a peak of the 
To — > S* channel. For r_, the dip and peak occurs in the 
T- ^ S and T_ T+ channels, respectively. Note that 
if the in-plane magnetic field is perpendicular to the dot 
main axis d (lower panel) , the relaxation channels for T_ 
and Tq do not vary at all, as the S — Tj^ anticrossing gap 
vanishes, AE = 0, and the exact compensation is trivial. 



Hyperfine versus spin-orbit induced relaxation 

Comparing the value of the nuclear and spin-orbit 
effective fields, we estimate the relaxation due to the 
former is typically three orders of magnitude smaller, 
{Bnuc/Bso)^ ~ 10~^, if the external field is of the order 
of Tesla. However, in a weakly coupled detuned double 
dot the nuclear spins can dominate over the spin-orbit in- 
duced relaxation in some cases, when Eq. (14) is satisfied. 
We plot in Fig. [6] the spin relaxation rates enabled by 
spin-orbit and hyperfine coupling, respectively. Panel a) 
gives the relaxation rate of the first excited state. The hy- 
perfine coupling becomes relevant only close to the S—T^ 
anticrossing along the easy passage. Here, the wide dip 
is narrowed (red versus the blue curve). However, the 
rate remains reasonably low, such that the easy passage 
survives. Adding the nuclear dominated area to Fig. 4 
a) of the main text would barely be visible. Panel b) 
shows the rate of Tq. We find that the hyperfine- induced 
relaxation is dominant for any in-plane magnetic field 
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FIG. 6. Calculated spin-orbit induced relaxation rates for 
an in-plane magnetic field orientation parallel (black curves) 
and perpendicular (blue curves) to the dot main axis d. The 
red curves show the hyperfine-induced spin relaxation, (a)- 
(c) Weakly coupled double dot (T = f ^eV) as a function 
of detuning for B = 2 T. The panels display the relaxation 
rates for the first excited state, the unpolarized triplet, and 
T- respectively, (d) Unbiased double dot as a function of 
interdot distance (in units of lo) for B = 5 T. The relaxation 
rate of To is shown. 

orientation if the unpolarized triplet is close in energy 
to the first excited singlet, as shown in Fig. 4 b) of the 
main text. Panel c) displays the relaxation of T_. At 
the S — T_ anticrossing, the spin-orbit induced relax- 
ation strongly peaks unless the in-plane magnetic field 
orientation is perpendicular to the dot main axis. At the 
anticrossing, also the hyperfine-induced rate is enhanced. 
Displacing the magnetic field from the easy passage, the 
spin-orbit rate quickly gains on magnitude, therefore the 
nuclear-dominated area on Fig. 4 c) would cover only 
a single point at its current resolution. In panel d) we 
show the relaxation rate for an unbiased dot. We choose 
To as an example, the state which is most prone to have 
relaxation dominated by nuclear spins in the biased dot. 
Here, the relaxation due to the spin-orbit coupling is sev- 
eral orders of magnitude larger than due to the nuclei, 
for any orientation of the external field. We observe a 
similar difference in rates for other states in this setup as 
well. 



Easy passage exploitation examples 
Dynamical nuclear spin polarization 

We sketch two schemes of dynamical nuclear spin 
pumping in Fig. [T] The first is the one originally pro- 
posed by Reilly et. al. in Ref. [35]. Here, the double dot 




FIG. 7. Sketches of two schemes of dynamical nuclear spin 
pumping. The blue arrows indicate the path which the system 
state follows during one pump cycle. 



is initialized in the 5'(o,2) state. Then the system is adi- 
abatically brought through the anticrossing (step 1), by 
which a nuclear spin is flipped, assuming the anticross- 
ing is due to the (transverse component of) the nuclear 
effective field (and not due to the spin-orbit coupling). 
Placing the system into the easy passage, which was not 
done in the experiment, thus offers improved scheme effi- 
ciency. The cycle is finished by resetting the system into 
the >S'(o,2)i achieved by a fast transition (step 2) and a 
subsequent relaxation (step 3). 

We propose here a non-adiabatic version of the scheme, 
which does not require the increasingly slower gap 
crossovers, by ending the step 1 inside the anticrossing. 
This necessarily requires to monitor the anticrossing po- 
sition, which however is possible. Again, the scheme is 
most efficient if the spin-orbit contribution to the an- 
ticrossing gap is minimized, what happens in the easy 
passage configuration. 



Spin polarization detection 

Here we propose a device which allows to detect the 
spin polarization of a lead using a weakly coupled dou- 
ble dot. The dot is connected to source and drain leads 
such that the current passes only through the left dot. 
The system is biased such that only (1,1) and (0,2) oc- 
cupations are allowed if the right dot is brought below 
in energy and the current is allowed to flow. The sys- 
tem is periodically brought above both source and drain 
leads so that the right dot is emptied. After such a reset, 
when the right dot is lowered in energy by a gate, the 
electron which is traversing the device may tunnel into 
the right dot and becomes trapped (we assume its spin is 
preserved). Now, if an electron with the same spin orien- 
tation enters the left dot, a T_|_ (i i) triplet state is formed, 
with a long lifetime such that the electron in the left dot 
tunnels out and the current flows. If, on the other hand, 
a spin opposite electron enters the left dot, the system 
quickly collapses into the 5(0,2) state and the current is 
blocked until the reset (see the sketch in Fig. [8]). As the 
result, the higher the spin polarization of the electrons in 
the source lead, the higher the current on average. The 



8 



J© © 

"'+.(1.1) 



b) 




FIG. 8. Spin polarization detection scheme, (a) The double 
dot is in a T.^ state and the current is enabled, (b) For 
the 5(0,2) state the current is blocked. 



scheme requires 



rTo-s-S ^ rioft-).lcad ^ Tt^-j-S, 



(25) 



where F are the rates for transitions corresponding to the 
indexes. The conditions in Eq. (j25p are only achievable 
in the easy passage configuration while r^^^s must be 
dominated by the hyperfine-induced relaxation. 
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